Measurement of the branching fractions and mass spectra for r lepton decays 

including Kg at Belle * 
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We report a study of r lepton decays involving if^ with a 669 fb~^ data sample accumulated with 
the Belle detector at the KEKB asymmetric-energy e^e~ collider. The branching fractions have been 
measured for the r" ^ -K-Rlur, K-R^iy^, n-K'^n°u^, K-K^Tr°Vr, -k' K%K°sVr and tt" A'§7i'g7r°z/^ 
decays. We also provide the unfolded mass spectra for t~ — >■ ir~ K%-k^Vt and r~ — >■ K~ Kgiv'^UT. 

PACS numbers: 



INTRODUCTION 



r lepton is the only known lepton massive enough to 
decay into hadrons. Its hadronic decays are ideally suited 
to investigate the hadronic weak currents. Former stud- 
ies for semileptonic r decays done at LEP and CLEO [1,] 
provide accurate results for measurements of branching 
fractions as well as spectral functions for one- or three- 
prong decay modes with any number of neutral mesons 
that do not suffer from Cabibbo and/or phase space sup- 
pression. However, because of limited statistics, their 
studies for Cabibbo-suppressed or multi-prong r decay 
modes do not provide sufficient information for investi- 
gating the hadronic structure or for testing important 
standard model parameters, e.g., \Vus\- 

In this analysis we use a data sample of 669 fb~^ cor- 
responding to 616x10^ T~^T~ pair events, which is two 
orders of magnitude larger than those that were avail- 
able prior to the S-factory experiments. The r+r~ pair 
events, accumulated with the Belle detector are pro- 
duced at the KEKB asymmetric-energy e+e" collider Q 
running in the energy range of the T (45") with a cross sec- 
tion comparable to that of BB [ij . It is also worth noting 
that the r+r^ produced at the T(4iS') are in a very clean 
experimental environment with low background. 

The leptonic decay modes are used here for tagging r 
events for a precise measurement of the four decay modes: 



and r K Kgir^Vr- The other two decay modes, 



and 



however, 



are selected using one-prong decays for tagging to in- 
crease the signal. 
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EVENT SELECTION 

The selection of t^t^ events is performed in two 
stages. At the first stage, loose conditions are applied 
for the primary r+r^ selection. A r+r" event produced 
back-to-back in the e^e~ center-of-mass (CM) frame is 
divided into two hemispheres according to the plane per- 
pendicular to the thrust axis. Each event is required to 
have only one track on one side of the hemisphere (tag 
side) and any number of tracks on the other side (signal 
side). 

As a second step, we classified the events into two 
types, leptonic and hadronic events, using the combina- 
tion of the number of charged tracks and their flavors on 
each side. The events on both sides containing a single 
track are classified as candidates of leptonic events, while 
the events where the tag side contains a single track with 
three or more tracks on the signal side are classified as 
hadronic events. 

Particle identification (PID) is used to determine the 
flavor of charged tracks on both sides by defining and 
using the likelihood variables to identify the electron and 
muon track on the tag side and the charged pion or kaon 
on the signal side except the daughter tracks. 

The leptonic events in which the t lepton on one side 
decays to e~i>(.VT and the other side to ii~i>^Vr (e-/i 
events), are selected for normalization. The branching 
fractions for t~ — ^ e^Pgi^r and t~ — fi~i>^iy-r are mea- 
sured with good precision and can be used as a normal- 
ization for the branching fraction measurement. The de- 
tection efficiency and its statistical error for e-fi events is 
(19.31 ± 0.03)%. 

Hadronic events are selected by different tagging with 
at least one Kg — >■ tt+tt", which is reconstructed with a 
vertex fit using the momentum of two oppositely charged 
tracks on the signal side. The lepton-tagged events are 
defined as those with one r decaying to leptons, while 
the other decays to hadrons. Lepton tagging is used for 
TT-ifg, K~K^, TT-Kl-K^ and K-K^gir^ events. Since the 
signal yield for events involving more than one is 



2 



smaller than that involving a single K^, we use both 
lepton and hadron tagging for TT~K'gK'g and ■k~ K'gK^-K^ 
events. The decay modes used for hadron tagging are 
T~ h~{>0 tt'^)i't where h = Tr,K. 

The signal ir^ 77 is reconstructed from the invari- 
ant mass determined from the momenta of two photons 
detected on the signal side. The distribution of the dif- 
ference between the invariant mass of the two photons 
and the nominal tt'^ mass normalized to the resolution, 
— m^o)/cr-y-y, is used to determine the number 
of genuine tt^'s and to estimate the level of background 
from mass sidebands. 
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TABLE I: Efficiency migration matrix £ij for six decay chan- 
nels. The first four rows show the efficiency for lepton tagging, 
while the last two rows show the one for lepton and hadron 
tagging. 



DETERMINATION OF BRANCHING 
FRACTIONS 

The branching fractions for decay modes containing a 
single K'g are obtained using 



B{t- 



N 



Sig 
X-l 



Ne-^ Be 



B„ 



(1) 



rSig 



where X is the signal channel under study, N^_i — 
Nx-i/^x-i is the number of events with an efficiency cor- 
rection, where the r lepton on one side decays into a 
signal mode and the r lepton on the other decays into 
leptons, Ne-fj. = N°^^/te-^j. is the number of events with 
an efficiency correction for e-/i events. Be and are the 
world-average branching fractions for — >■ e^DeVr and 
Q, respectively. 
In order to increase the number of candidates, 
■K~ K^K^ and tt' K^K^n^ events are selected by both 
lepton and hadron tagging. To determine the branching 
fraction for multiple K'g modes, we used 



B{t- 



N' 



Sig 
X-A 



2 X Nrr X B{t- A-VV) ' 



(2) 



where A indicates a one-prong combination of particles: 
e.g., — > l^viVr and — > /i^(> vr*');^,-. 
In both cases, the number of signal events N,^ is given 



Sig _ 



Data 



(3) 



where i is the true decay channel of interest and j is the 
decay channel reconstructed in the analysis. Nf"^^^ is the 
number of selected events and N^^ is the estimated back- 
ground coming from decay modes other than six under 
consideration, and non-r processes. £ij is the efficiency 
matrix. The diagonal components are the efficiency for 
the i-th channel and the other components are the prob- 
ability of migration from the i-th to j-th decay modes. 

For the branching fraction measurement, the number 
of the background events coming from r decays other 



than the six modes studied should be determined. To 
determine the background, generic MC samples of r lep- 
ton decays generated using branching fractions provided 
by the Particle Data Group (PDG) are used. The non-r 
decay contributions are relatively small, at a level of 1%, 
and dominated by qq continuum events. 



EFFICIENCY AND CORRECTIONS 

In order to determine the signal efficiency and cross- 
feed rates for the decay modes of interest, a 17.2 million 
T+r^ MC sample is used. Since a difference between MC 
and data is inevitable, we evaluate corrections associated 
with a difference between MC and data for FID and the 
reconstruction of the neutral particles, and tt". The 
corrections of the FID efficiency for electron and muon 
are derived from two-photon events and those for charged 
pions and kaons are derived from D*~ — > D^Tr~ and 
— > K^TT^ control samples. The correction for the Kg 
reconstruction efficiency obtained is 0.979 ± 0.007 using 
D*- D°Tr-, D" K^TT+TT- control samples. The 
MC-data correction for the 7r° efficiency is determined to 
be 0.957 ±0.015. 

For the tt^ Kgir'^ and Kgir'^ modes, the number 
of cross-feed events from tt^K^ and K^K^ with spu- 
rious TT^ is estimated and subtracted by using sidebands 
of S-yj. By this subtraction, the cross- feed contributions 



from TT^ KgUr and/or — > K^Vr are removed 

successfully. 

By taking into account the corrections and the effi- 
ciency losses by subtraction, the efficiency migration ma- 
trix £ij is obtained and summarized in Table HI The di- 
agonal elements represent the signal efficiency, otherwise 
- the cross-feed rates among the decay modes of interest. 



SYSTEMATIC UNCERTAINTY 

The sources of systematic uncertainties are catego- 
rized into the uncertainty of detection/reconstruction ef- 
ficiency (EFF), hadron decay models (HDM), background 
estimation (BGE), normalization (NORM) and rejection 
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AB/B {%) 

source -n'Kl K- -k- K^-k'^ K-K'^tt° vr-^gifg tt' K'^.K^.tt'^ 



EFF 2.7 3.5 3.2 3.3 3.4 5.6 

HDM - - 0.3 3.4 

BGE 0.2 0.3 1.9 0.4 1.8 3.2 

NORM 0.4 0.5 0.4 0.5 1.3 1.3 

7 veto 0.1 1.8 1.2 1.5 1.0 2.0 

Total 2.8 3.9 4.0 5.0 4.2 6.9 



TABLE II: Summary of the systematic errors. 

iT^Kg i -0.28 -0.08 OOl -0.01 3.9c-3 

K~ K*^ 1 0.03 -0.14 -0.9C-3 0.3c-3 

tv-K'^-k" 1 -0.21 -0.06 0.02 

K~K^Tr° 1 2.9e-3 -1.0o-3 

1 -0.37 

77-ifgK°7r° 1 



The correlations of these uncertainties are also taken 
into account by the covariance matrix cov{Bi,Bj) us- 
ing the formula given in Ref. The results 
including both systematic and statistical uncertain- 
ties are shown as the coefficients of the covariance 
cov(S, ,Bj)/ y^cov{B^,B,) coviBj,Bj) in Table Hill 



RESULTS 

Using the formula given in Eq. ([T]), Eq. ^ and 
Eq. ([3]), the branching fractions for six decay modes 
are obtained simultaneously. The branching fraction for 
T- n-Klur is found to be (4.13 ± 0.01 ± 0.12) x lO^^^ 
which is consistent with our previous study The 
branching fractions for K^Kgh'r, tt~ Kgir'^ i/r and 

K-K^Tr°Vr are measured to be (7.36±0.04±0.29) x lO""^, 
(1.92±0.02±0.08) X 10-3, and (7.44±0.11±0.37) x lO^*, 
respectively, which are among the most precise results 
to date. The branching fractions for n'KgKg and 
TT~ KgKgTr'^ have been measured for the first time to be 
(2.39±0.03±0.09) x lO"'^ and (2.06±0.13±0.14) x lO'^, 
respectively. 



STUDY OF T- TT-K^sK'^n^^iyr 

Not only the branching fraction but also the hadron 
current can be studied in t~ — > '!r~ K^Kgir^Vr using the 
invariant mass oi KgKgir'^ , M{KgKgTT^). As can be seen 
in Fig. [U the M{KgKg7T'^) distribution shows a signifi- 
cant peak at around 1280 MeV/c^, which can be under- 
stood as the /i(1285) resonance in r decays. In order to 
determine the mass and width for the resonances, we used 
the probability distribution function (PDF) consisting of 
a Breit-Wigner function for the signal and a second-order 
polynomial for the background. So the signal is parame- 
terized as 

a( ^ ^ mMTp 

(M2 -m2)2 + {mToy ^ ' 

where M and Fq are the nominal mass and width, m is 
the mass of the ifgif^Tr^ system. The total number of 
signal events is found to be 80±16. The mass and width 
are 1278±4 MeV/c^ and 35±9 MeV/c^, consistent with 
the /i(1285) properties. From the total number of — ^ 
Tx-KlKl-K^Vr and the number of /i(1285) K%K%Tfi 
among the — > KgKgir'^h'T events, the branching 
fraction of — ^ /i(1285)7r^i/^ — > K'^K^TT^i/r is 
found to be (1.05 ± 0.24) x 10"^. In this fit, possible 
interference between /i and non-/i amplitude is ignored, 
which will cause considerable model uncertainties in the 
/i parameters and the branching fraction. Further stud- 
ies are in progress. 



TABLE III: Coefficients of the covariance matrix for statisti- 
cal ® systematic uncertainty covariance. 

of energetic photons (7 veto). The uncertainty of effi- 
ciency consists of several items: track finding, FID, Kg 
and TT^ reconstruction and the ir'^ sideband subtraction 
for TT^Kgir'^ and K^Kgir'^. These uncertainties are inte- 
grated into the efficiency migration matrix jsj. The un- 
certainties from the decay model implemented into the 
T MC sample are estimated by using weighting sets of 
MC samples generated with different hadron decay mod- 
els. The background uncertainty mostly contributed by 
other r lepton decays is estimated by varying the frac- 
tions according to their uncertainties in PDG. By adding 
these uncertainties in quadrature, the total systematic 
uncertainties for tt-R^, K^K%, tt~ KItt° , K' KItt^ , 
TT-K%K% and tt' KIK%tt^ are 2.8%, 3.9%, 4.0%, 5.0%, 
4.2% and 6.9%, respectively. The summary of systematic 
uncertainties is shown in Table HIl 




FIG. 1: Invariant mass of K%Kln° in t- ^ KIkItt°v^ 
events. A significant /i(1285) signal is seen. The solid hne 
is a fit with a Breit-Wigner and a second-order polynomial 
function. 
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FIG. 2: Unfolded invariant mass spectra for tv~ K^tt^ (Mt) in TT'K"gTv°Ur and AT^J^gyr" (right) in K'Kl-K°Vr. The 

error bars show statistical error and are within the closed circles in many points. 



UNFOLDED MASS SPECTRA FOR 

T- TV-Kl-K^Vr AND T- ^ K^-K^ Vr 

The mass spectra provide important information for 
hadronic weak currents: the tt^ K^i:^ mode contributes 
to the Cabibbo-suppressed weak current, while the 
K^Kgir'-' mode has both vector and axial-vector com- 
ponents. The axial-vector component is interesting as 
the Wess-Zumino-Witten anomalous term [91] . The mea- 
sured mass spectra are distorted due to the finite res- 
olution and the efficiency of the detectors. In addition, 
there is a sizable cross-feed between these two modes. By 
using the singular value decomposition (SVD) unfolding 
technique, the total mass and submass spectra are ob- 
tained from T~ — > TT^ Kgir'^ and K^tt^Vt by 
taking into account the effects from the detector efficien- 
cies, resolutions as well as the cross-feed backgrounds. 
The results are shown in Figs. [D A clear resonance-like 
structure due to _ft'i(1280) is seen in — > tt" K^tt^Vt. 



CONCLUSION 

Using 616 X 10^ T+r~ events coUlected with the 
Belle detector, we measured several branching fractions 
for hadronic t decays involving K^: TT~Kg, K~Kg, 
TT-if^Tr", K-KItt°, T^'K^Kl and tt- KIkItt^ and pro- 
vided the covariance matrix, which includes both sta- 
tistical and systematic uncertainties. We also provide 
the correlation of the systematic errors as the covari- 
ance matrix. The invariant mass of K'gK^i:^ in the 
'K~ K^K'gir'^ mode shows a clear peak around the mass 
of 1280 MeV/c^, which most probably comes from the 
/i(1285) current. In addition, we measured the unfolded 
mass spectra for r~ tt' K^n^Vr andr" K^ir^Vr 
using the SVD unfolding technique. 
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